Effective action approach to the Leggett's mode in two-band superconductors 
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We investigate a collective excitation (Leggett's mode) corresponding to small fluctuations of the 
relative phase of two condensates in two-band superconductor using the effective "phase only" action. 
We consider the possibility of observing Leggett's mode in MgE>2 superconductor and conclude that 
for the known at present values of the two-band model parameters for MgE>2 Leggett's mode arises 
above the two-particle threshold. 
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I. INTRODUCTION 

The problem of collective modes in superconductors is almost as old as the microscopic theory of superconductivity. 
Bogolyubov Q and Anderson ]!| discovered that density oscillations can couple to oscillations of the phase of the 
superconducting order parameter via the pairing interaction. In a neutral system these collective soundlike oscillations 
are called Bogolyubov- Anderson- Goldstone (BAG) mode. In a charged system the frequency of this mode is pushed 
up to the plasma frequency due to the long-range Coulomb interaction |g| . 

Despite the fact that the basic physics of the BAG and plasma modes in superconductors is well understood many 
years ago, it appears that a more modern approach provides even better insight to the origin and universal character 
of these modes allowing also to tackle less settled questions. A main idea beyond this approach is rather simple: since 
the collective modes present usually low energy degrees of freedom, to study them it is sufficient to have only an 
effective Lagrangian (or action) that describes low-frequency, low-wavelength dynamics of the phase 9 of the pairing 
field instead of working with the original Hamiltonian of the system. (Note that although plasmons are in general 
high-energy excitations, they still can be treated using an effective theory 0j.) 

The most convenient way of deriving such a Lagrangian is to change the variables (see, for example, W) for the 
complex pairing field $ — > A exp(i0) with real modulus A and phase 9. The simultaneous transformation for the fermi 
field ipo- — > \o- exp(i9/2) allows from the beginning to separate in the Hamiltonian the only phase degree of freedom 
relevant for the effective theory. In the theory of superconductivity this replacement of the variables that resembles 
a gauge transformation has, in fact, a long history and it was probably firstly used in Q and then in || (see also 
Rcfs. in the review 0). More recently this transformation was used by many authors. For example, in g it was 
used to study the problem of the Galilean invariance of the effective Lagrangian at T = 0. The finite temperature 
time-dependent effective actions for the phase field in s-wave Q and d-wave |l0| neutral superconductors were derived 
addressing an old problem of time-dependent generalization of the GL theory. It is also convenient to investigate the 
plasma mode within this formalism taking into account the long-distance Coulomb repulsion between electrons |l|, |^] 
(see also more recent papers |y], [T^] ) . 

Besides these rather commonly known modes there are also other collective modes such as Carls on- Goldman mode 



|L3| which can appear when the Coulomb interaction gets screened and Leggett's mode 14 which is present in two- 
band superconductors. Recently Carlson-Goldman mode has been studied in d-wave superconductors Q and colour 
superconductivity |ll| using the approach based on the modulus-phase variables, but to the best of our knowledge 
Leggett's mode has not been considered yet using this transformation. Thus it would be interesting and instructive 
to apply this method to obtain Leggett's mode. 

Physically Leggett's mode [|l4| is a collective excitation corresponding to small fluctuations of the relative phase of 
two condensates in a clean two-band superconductor. There are some indications that MgB2 superconductor with 
T c = 39K [jljj discovered about one year ago is an example of such a two-band system. Thus the main results of the 
present paper are two- fold and can be summarized as follows. 
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i) Leggett's mode is obtained using the modulus-phase variables in the path integral formalism; 

ii) It is considered whether this mode can be observed in MgB 2 superconductor. 

The paper is organized as follows. In Sec. |l| we introduce the two-band model and represent the partition function 
of the system using the Hubbard-Stratonovich transformations for pairing and Coulomb parts of the interactions. In 
Sec. [II we express the effective thermodynamical potential of the system in the modulus-phase variables and obtain 
the system of equations for the superconducting gaps. Then we separate the part which describes the collective phase 
modes. In Sees. |iy| and ^ we derive the dispersion law for Leggett's mode in neutral and charged systems, respecti vely. 
Sec. VI we investigate the possibility of observing Leggett's mode in MgB2 superconductor. We conclude in Sec. VII 
with a discussion and summary of our results. 



II. MODEL TWO-BAND HAMILTONIAN AND HUBBARD-STRATONOVICH TRANSFORMATIONS 



Let us consider the following action (in our notations the functional integral is expressed via e ) 
rP r 2 



S = - I dr E E / d2r V>l(T, r)0rlMT, r) + H(t) 
Jo ,_, „ J 



where the Hamiltonian H(t) is 

2 



T 



(1) 



H( T ) = EE / d 2 ripUr,r)[ Si (-iV) - fx]^(r,r) 

1=1 a J 
2 

~ \ E E/ dri / dr 2V'i CT (r,r2)^(T,ri)F i j(ri;r 2 )^(r,r 1 )^v(r,r 2 ) ^) 

+ \ j dYl J dT2 ( E V't(T,ri)V'i<r(T,ri) -n J 14(r x - r 2 ) I ^ ^t CT ,( r ,r 2 )^ CT '(T,r 2 ) - 

Here ipi<r(T, r) is a fermion field with the spin a =1", J,, a = — a, i,j = 1,2 is the band index, £i(k) = k 2 /2mj is the 
dispersion law in i-th band with the effective mass of the carriers m;, r is the imaginary time and Vy (ri; r 2 ) = Vij5(ri — 
r 2 ) is an attractive short-range potential, V c (yi — r 2 ) is the long range Coulomb interaction, n is the neutralizing 
background charge density. Throughout the paper we call the superconducting system neutral if the last term of 
Eq. (|2j) is omitted and charged if this term is taken into account. Even in the latter case the whole superconductor 
remains, of course, neutral due to the neutralizing ionic background. Throughout the paper h = ks = 1 units are 
chosen. 

Introducing Nambu variables 



^(r,r) 

we rewrite the action as a sum 



S — Sq + Sp a i r + Sc 



of free 



pairing 



So 



r p n 2 
Jo J i=1 

fl 3 f 2 
Spair = J drjdr^ V lJ * l t T + * l (a;)*]r_* i (a;), 



and Coulomb 



Sc = -\f Q dT J dr i / rfr 2 ^X>fc)7S*i(aO-nJ V c (n -r 2 ) (^^(^3%^) 



(3) 



(4) 



(5) 



(6) 



(7) 



3 



parts. Here r± = (tj ± ir%)l% t\ (A = 1, 2, 3) are Pauli matrices. 

The easiest way to treat S pa i r is to introduce Hubbard-Stratonovich fields $i for each band, so that 

Spuria,*;, *t) = J^dr J dr [-su|$i(x)| 2 - 0, 2 |*2(x)| 2 + ffl2 (^0r)$ 2 (z) + $5(z)*i(s)) 

2 



i=i 



where the coupling constants are expressed in terms of the original constants Vij 



9ii =V n [l 



V 2 

v 12 



922 = V 2 2 1 - 



12 I 1/ . h 12 



, 512 = V12 1 - 



^11^22 



V\\V 22 J V V11V22 

For the repulsion part only one Hubbard-Stratonovich field ip is necessary: 

= J dr J dri J dr 2 ~e(p(r, ri)V c _1 (ri - r 2 )e(p(r, r 2 ) + ieip(r, n) *J(a;)T 3 *i(x) - <5(n - r 2 ) 
Thus the partition function of the system can be presented as 



(8) 



(9) 



(10) 



(11) 



III. EFFECTIVE POTENTIAL IN THE MODULUS-PHASE VARIABLES 

A. Modulus-phase variables 

The modulus-phase variables Aj and 6i in two bands are introduced exactly as discussed in |7j] : 

sl 9 i (r,r)/2 q 



$i(r,r) = Ai^^exp^i^r)], *i(r,r) 







-i9i(T,r)/'< 



r,r)/2 Ti(T,r), 



(12) 



making the moduli A, of the pairing field real. We note that since we restricted our analysis to T < T c in zero 
external field, so that no vortices are present, the phases 9i(r, r) are non-singular, so that the transformations ( |l2] ) 
do not change the magnetic field "seen" by quasiparticles. Now absorbing gn and 522 in Ai and A 2 and integrating 
out neutral fermions Y.; we obtain (see e.g. M, [l5| [l9| ) 



Z = J A j; PA 4 m^exph/3^(A 4 ,0 4 ,^)] 



(13) 



with the effective thermodynamical potential 



pCl(&iA,<P) = / drdv 
Jo 



A 2 , Ai 2g 



12 



drdridr 2 



9ll 922 911922 

1 



AiA 2 cos(6>i - 2 ) 



e<^(T,ri)V c 1 (ri,r 2 )e< y 9(T,r 2 ) + ieip(r, ri)n<S(ri - r 2 ) 



Tr LnGr 1 -Tr LnG^ 1 , 



where the Green's function 



G7 1 = - MdOi) = -Idr + r 3 ( — + M ) + nAifr, r) - Mdd t ) 



with 



idrBi . . . J (V0 4 ) 2 



- / 



iV 2 ^ iV0i(r,r)V 



4to,- 



2to ? - 



(14) 



(15) 



(16) 
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that depends only on the time and space derivatives of 9i, but not on the phase 9i itself. Then we can represent f2 as 
the sum 



where 



fi(Ai, 6i, ip) ~ fi k in(Ai, d9 t , ip) + n po t(A i , 6i, (p) 



2 oo 



i=l 71=1 

is the sum of the energies of the phase fluctuations in each band and 



(17) 



(18) 



Ai— const 



fi P ot(Ai, 9i, if) = 



A? , A| 



311 322 311322 



2 ' /12 - Ai A 2 cos(0i - <9 2 )] - Tr Lngf 1 - Tr 1 



Ai— const 







d,Td,YidY2 



-eip(T,r 1 )V c 1 (ri,r 2 )e< y £>(T,r 2 ) + ie(^(r, ri)n^(ri - r 2 ) 



(19) 



In f2p t the most important for the appearance of Leggett's mode term is the Josephson coupling energy of the 
condensates in two bands. This term explicitly depends on the relative 6\ — 9 2 phase of two condensates. 



B. The system of gap equations 

As pointed out in |TJ] while V\\ and V 22 are completely fixed by the physics of the problem, the phase of V\2 is a 
matter of convention (corresponding to the choice of the relative phases of the Bloch waves in the two bands). As in 
we choose Vvz to be real and positive. In this case the condition of minima of £l po t with respect to 9\ — 6> 2 gives 
9\ = 92, so that the system of the gap equations dtt po t/dAi = has the form 



JAx - §^A 2 - A l5ll iV 1 F(A 1 ) = 0, 
|A 2 - f^Ax - A 2522 iV 2 F(A 2 ) = 0, 

where JV» = rriipFi/ (2tt 2 ) is the density of states (per spin) in z-th band (pfi is the Fermi momentum) and 



(20) 



F(AQ=r , * tanhV^ (21) 
V t} Jo VeTAf 2T v > 

with the Debye frequency, lojj which is for simplicity assumed to be the same in each band. 

This system of the equation can be transformed to the standard form derived by Moskalenko pp| and Suhl et al. 

f Ai[l - VuNiF(Ai)] = A 2 V 12 N 2 F(A 2 ), 

1 A 2 [l-y 22 A 2J F(A 2 )] = A 1 ^ 12 A 1 F(A 1 ). 1 ' 



It appears, however, that in contrast to Leggett's approach |L4j there is no need of explicit use of the system (|22j) to 
obtain the spectrum of collective excitations. In what follows we assume that the values of the superconducting gaps 
Ai(T) and A 2 (T) can, in principle, be determined from the system (123) or taken directly from the experiment. 



C. Effective potential for the collective modes 

Since the values of the gaps are fixed, i.e. there are no amplitude fluctuations that is reasonable for T <C T c it is 
sufficient to consider only the part of effective potential depending on 9i and ip. 

Following the same route as described, for example, in Refs. ]9|, [l0| [l{|, we arrive at the following frequency- 
momentum representation 

8 n^ooJ ( 27r > V 311322 
2 

+ e<p(-K)]Mr l 



(23) 



6i{K) 
e<p(K) 
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where the Josephson term was expanded up to quadratic term and unimportant constant was dropped out. The 
matrix M.i in (23) is 



AC 1 = 



-2m n m 33 (K) + 2K a m^ {K) -w 33 {k) 



(24) 



and we introduced short-hand notations K — (ifi n ,K) with Q n = 2irnT and K being 3D vector (summation over 
dummy indices a, (3 = 1,2,3 is implied). In Eq. (|24|) l K a!3 = 4 Aq^ + "IIqo is the bare (unrenormalized by the phase 
fluctuations) superfluid stiffness, where the current-current polarization function, l IIoo, is 

°° r d 3 k 

i n^(ir!„,K)=T J j^**mWn,K;iw h \L)v Fia (k)v Fif} (lL) (25) 



i=-c 



with the Fermi velocity VpiaQs) = d£i(lt) / dk a \k=k Fi (here £i(k) = £j(k)— /x); the density-density polarization function, 

i n 33 , is 

°° f d 3 k 

TI 33 ^„,K)=T^ J ^—Vas^.Kjiwj.k) (26) 



i 

and the density-current polarization function, 1 TIq 3 is 

d 3 k 



(k). (27) 



V Ak in Eqs. (|25|) - (|27]) is given by 

WCifin, K; iw!,k) = tr[*e(»wi + *fi n , k + K/2)r x i a(^,k - K/2)t„] , (r = J) , (28) 
where the neutral fermion Green's function coincides with the usual Green's function of the BCS theory 

ire i,i iw n J + r 3 ^(k) -nAi 

g(iw n , k) = 2 T &n7\ _i_ A2 — ' w ™ = ^(2" + 1 ) T - ( 29 ) 

l Ag^ above is the first order contribution to the superfluid stiffness: 

V = I -^WkVr 1 ^ = ^<5 Q/3 , (30) 
J \2^) m t 



with 



^(k) = l-|^tanh^, B i (k) = v /C?(k)+Af. (31) 

Writing Eq. fl23] ) we omitted the time derivative term linear in the phase (see e.g. |l(| pJj ) which is irrelevant for the 
present analysis. 

For the purpose of completeness we give the explicit expressions for the polarizations (^5|) - (|2^) in Appendix [A|. 
These expressions are important if, for example, one is interested in the damping of the collective modes. In what 
follows we consider these polarizations in the hydrodynamic limit, f2„ = and K — » at T = 0. In this case 
TI 33 (0, 0) = -2Ni, l A a/3 (0, 0) = = n l /m t = 2N l v 2 Fl /3 and l U% 3 (0, 0) = 0. Recall that N t j s the density of states 
per spin, while in [fbil the density of states per particle was used. 

As was mentioned above, calculating the values *II 33 (0,0) and l A a ^(0, 0) there is no need to substitute the gap 
equations ( p2|) since the coupling constants Vij do not enter z II 33 (0, 0) and l A Q,3 (0, 0), but enter only in the Josephson 
coupling term. By contrast, the original derivation of Leggett |l4[] (see after Eq. (3.8)) explicitly uses the system 
Note that in a more simple case of one band system the Josephson term is absent and one immediately gets the BAG 
mode without referring to the gap equation || . 

We note also that while ^ ^m is zero for the case considered in the present paper, this term is crucial for the existence 
of Carlson-Goldman mode pl| |l5| . 

Having the general representation for the thermodynamical potential ip} we are ready to obtain the spectrum 
of collective excitations. We start with a more simple case of a neutral superconductor. 
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IV. COLLECTIVE EXCITATIONS IN A NEUTRAL SUPERCONDUCTOR 



To consider neutral superconductor one can formally set e = in (23), so that the terms with the electric potential 
ip disappear from the equations and we arrive at 



8Vi 



12 



) 3 \VnV22 - 
J29 l (-K)M ei 1 9 t (K)) , 



A 1 A 2 (9 1 {-K) - 9 2 (-K))(9 1 (K) - 9 2 {K)) 



(32) 



i=i 



where 

m^ 1 = -nlm^R) + i A a,3 (K)K a K p - tn n K a m« 3 (K) - m n Kjn% Q {K). 

As was already mentioned, we consider Mg i in hydrodynamic (f2 n = 0, K — ► 0) limit at T = 0: 

Mg i L = 2N i (n 2 n + c 2 K 2 ), 
where cf — is the velocity of the BAG mode in i-th band, so that Eq. ( |32"| ) becomes 

v{K) 



T 00 P d\i 



i{K) 



with 



e- 1 



+ A -A 
-A Ml ' 



A 



A 



8T/i 2 AiA 2 

V U V22-V&' 



(33) 
(34) 

(35) 
(36) 



Finally, solving the equation detO 1 = for collective modes and making an analytical continuation iCl n — > uj + iO 
we arrive at 



1 



with 



o 2 + (cj + c 2 2 )K 2 ± Jut + (c 2 4) 2 * 4 



2 _ Ni +N 2 8T/i 2 AiA 2 



Ni -N 2 

m + N 2 



2N X N 2 V n V 22 -V 2 2 



(37) 



(38) 



(For a direct comparison with [Q recall that the density of states used here is twice less.) In the limit K — > 
{vFiK <C wq) one obtains from (|37|) 



2 2 jv-2 z 

w = c K , c = 



lu' z = OJq + v 2 K 2 , v 



2 _ A^ie 2 + N 2 cl 
Ni+N 2 
N lC 2 + 7V 2 c 2 
Ni +N 2 



2 for 
for 



(39) 



The first solution of ( |39| ) corresponds to BAG mode, while the second solution is Leggett's mode. This collective 
mode is only possible if ui 2 > 0. Since V% 2 > this implies that Leggett's mode exists for Vi%V 22 — V 2 2 > 0. 



V. COLLECTIVE EXCITATIONS IN A CHARGED SUPERCONDUCTOR 



As was discussed in Introduction the long-distance Coulomb interaction has a drastic influence on the BAG mode 
transforming it in the plasma mode. Here we consider how Leggett's mode is affected by the Coulomb interaction. 
Let us rewrite action (|23|) in the hydrodynamic limit as one matrix 



T 



E 



n— — oo 



dK 

(2^)3 



[01 (-#) 9 2 (-K) e^i-K^M- 1 



9i{K) 
9 2 {K) 
eif(K) 



(40) 



7 



where 



AT 



Ma 



A 



-A 



-A 



Ml 



A 



UVL n N 1 UQ n N 2 4(2iVi 



-4ifl n Ni 
-4ifl n N 2 
2N 2 + V- 1 {K)), 



(41) 



and M e and A are given by Eqs. (34) and (p6|), respectively. As one can readily see, setting e = in Eq. 



reduce the action (|o|) with 3x3 matrix M.~ x that describes the phase fluctuations in the charged system to the 
action ([35]) with 2x2 matrix O -1 that describes the collective excitations in a neutral system. 

Again the spectrum of collective excitations can be found solving the equation det-M -1 = 0. Neglecting the inverse 
of Coulomb interaction in the limit K -»• (K _1 ( K ) = K 2 /(47re 2 ) -> 0), we obtain 



n 1 



(JVi + N 2 )c\cl 



N 2 cl 



(42) 



Since the inverse Coulomb potential ~ K 2 is neglected, the equation for the collective modes has the only solution 
describing Leggett's mode. Comparing ( |39| ) and ( [t2"| ) one can see that in accordance with Q the value of ujq is 
insensitive to the Coulomb interaction but the velocity v is affected. It is interesting that the value of v itself does 
not depend on the explicit form of V C (K). 

Including the Coulomb interaction V^(K) one can also check that the equation det.M _1 = has a plasma mode 
as a solution. Indeed considering the simplest V12 = {A = 0) case we obtain the plasma frequency for in-phase 
oscillations in both bands 



nt = 8ire 2 (N lC j 



N 2 c 2 2 ) 



Aire 



Hi 



«2 



mi m 2 



(43) 



VI. IMPLICATIONS FOR MgB 2 



There are currently many indications that the recently discovered MgB2 superconductor Jl7| can be described by 



the classical two-gap model [|20| 
In particular, in 24 



21 which convincingly fits the specific heat |22| and penetration depth measurements 
even the values of the coupling constants for the system of equations ( |22] ) for the 
two-gap and two-band model of MgB2 were given. 

We note that one of the alternative explanations of the observed anisotropy of the upper critical field uses a model 
with anisotropic s-wave pairing [ p7| which is also rather close to the two-gap and two-band scenario. 

To be experimentally observable Leggett's mode should have the value of luq in p2] ) well separated from the two- 
particle threshold given by the smallest gap, Ai. Here we estimate the value of ujq using recently suggested values 
of the coupling constants that enter the system of equations (p2|). Introducing the dimensionless coupling constants, 
Xij = NiVij that are often used for the description of the two-band model, we may rewrite Eq. (|38| ) in the following 
form 



A12 + A2 



A11A22 — A12A21 



4AiA 2 . 



Our estimates of luq are summarized in Table |. 



(44) 



TABLE I: The estimate of the frequency wo- 
are used. 



The values of the superconducting Ai = 1.8 meV and A2 = 6.8 meV from 



Reference 


An 


A22 


A12 


A21 


luq, meV 


^o/2Ai 


Liu, et al. |24| and Barabash 2C] 


0.96 


0.28 


0.16 


0.22 


8.9 


2.5 


Golubov, et al. p[ 


1.017 


0.448 


0.213 


0.155 


6.5 


1.8 



They show that for the values of the two-band model parameters known at present for the two-band model of 
MgB 2 , Leggett's mode arises above the two-particle threshold, u /2Ai > 1, and unlikely to be observed. We do not 
exclude, however, that Leggett's mode can be observed in MgB 2 if the values of the coupling constants are revised, 
so that the interband coupling constants A12 and A 2 i would become smaller allowing lu /2Ai < 1. The observation 
of Leggett's mode would provide an additional insight to the underlying physics of this superconductor. 
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VII. CONCLUSION 



To conclude, we readdressed the collective excitations of the relative phase of the two condensates in a clean two- 
band superconductor using the effective "phase action" formalism. This formalism has proved itself as a convenient 
and economic method of the investigation of the collective modes in superconductors. Our estimates of the lowest 
frequency of Leggett's mode for MgB 2 show that it can hardly be observed in this superconductor. 



VIII. ACKNOWLEDGMENTS 



This work was supported by the research project 20-65045.01 by the SCOPES-project 7UKPJ062150.00/1 of the 
Swiss National Science Foundation. 



APPENDIX A: THE EXPRESSIONS FOR POLARIZATIONS II 



The expressions for the polarization functions IT^ (for simplicity we omit the band index i) are (see e.g. |], 

u 33 (m n ,K) 



and 



d 3 k 
1 + 



111 


0- 


3 12 








E+ 



E^E, 



1 



£L + iCL 



V±(k) ee 



E+- E_+ ifl n E+-E-- ifl n 
"+" 



+ E-- itt r 
n F {E_) 



[1 - n F (E_) - n F {E + )] 
n F (E + )}\v±(k), 



v Fa (k)v F p(k), 
1, 



n? 3 (ffi„,K) 



f d 3 k 




J (2^) 3 


{(" 


+ 


£- 


\2E + 


2£_ 



2E 4 



2E_ 



1 



1 



E + + E- + iVt r , 



E+ 



1 



1 



E+-E-+ ia n E+-E-- ifl n 



E_ - m r , 

[n F {E- ) 



[l-n F {E-)-n F {E + )\ 
n F (E + )] \ v Fa (k), 



(Al) 



(A2) 



where ^± = £(k ± K/2), E± = E(k ± K/2). One can also check that U%(iQ, n , K) = n^ 3 (iO n ,K). The second term 
m Eqs. Q and @ gives the contribution of the thermally excited quasiparticles (i.e. "normal" fluid component). 
This is the term responsible for the appearance of the Landau terms in the effective action fi flOl . 
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